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Abstract 

In this paper we suggest an implementation of Runge’s method for 
solving Diophantine equations satisfying Runge’s condition. In this 
implementation we avoid the use of Puiseux series and algebraic coef¬ 
ficients. 


1 Introduction 

Consider the Diophantine equation 

F(x,y) = 0 

in integers x, y and where F is a polynomial with integer coefficients. We 
shall assume that F is irreducible over C. The equation F = 0 then rep¬ 
resents a geometrically irreducible algebraic curve, which we denote by C. 
Denote by g the genus of C and the number of branches at oo by s. From 
a well-known theorem of Siegel [Si] it follows that if s + 2g — 2 > 0, then 
the number of integer solutions to F(x, y) = 0 is finite. Recently P.Corvaja 
and U.Zannier [CZ] gave a very surprising alternative proof of this fact 
using W.M.Schmidt’s subspace theorem. Unfortunately, neither proof of 
Siegel’s theorem gives an algorithm to actually solve the general equation 
F(x,y ) = 0. Only in very special cases this is possible. For example, 
A.Baker’s method of linear forms in logarithms allows one to solve equa¬ 
tions of the form y q = f(x) for any given q and any / € Z[x] having three or 
more distinct zeros (see for example [ST, Ch 6]). Yu.Bilu [Bi] studied nec¬ 
essary conditions for the applicability of Baker’s method and found several 
new instances where F(x,y) = 0 can be solved in principle. 

In this paper we take up an old paper of Runge [Ru, 1887] where equa¬ 
tions of a particular kind are solved. As introduction consider the equation 
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F(x, y) = 0. Let d be the total degree of F and denote the sum of all terms 
of total degree d in F by Fq. Suppose that Fq factors as a product of two 
non-constant relatively prime factors Fq = GqHq. This is called Runge’s 
condition. The branches at infinity of C either correspond to Go or to Ho. 
Runge’s idea was to construct a polynomial P(x,y) £ Z[x,y], non-constant 
on C. such that P(x, y) —> 0 as we move to infinity along one of the branches 
corresponding to Go- For sufficiently large x,y the integer points on these 
branches should satisfy P(x,y) = 0 because P assumes integral values at 
these points. We can then find them by elimination with F{x, y ) = 0. Simi¬ 
larly we deal with the branches corresponding to Hq. This idea is described 
in the introduction of Runge’s paper, so an algorithm to solve the equation 
is in principle there. In the present paper we shall turn Runge’s idea into 
an actual algorithm that can be fed to a computer. 

In [Sch] we find a generalisation of Runge’s idea if one considers weighted 
degrees. We shall present this generalisation in a slightly different language 
using Newton polygons. Let us write F = n fm,nX m y n where the sum¬ 
mation extends over a finite set of integer pairs. For each pair m, n with 
f m ,n / 0 we draw a rectangle with vertices (0,0) and (m. n ) in the plane. 
The Newton-polygon of F is the convex hull of these rectangles. We de¬ 
note it by Np. The edges of Np, not contained in the coordinate axes, 
are called the slopes of Np. To every slope E of Np we can associate the 
sum Fp(x,y) = Yl(m n)eE fm,nX m y n . The theorem of Runge and Schinzel, 
presented in a slightly different form, reads as follows. 

Theorem 1.1 (Runge, Schinzel) Let the notation be as above. Suppose 
that the Newton polygon of F has either two distinct slopes, or one slope E 
and Fp factors into two nonconstant, relatively prime polynomials in Z[x, y\. 
Then the equation F(x,y) = 0 has finitely many solutions. 

Using the proof of this theorem it is possible to give explicit upper bounds 
for size of the solutions (x,y). This is done in [HS] or [W]. However, the 
Tipper bounds are so large, even for small parameters, that using them for 
an exhaustive search on x, y is impossible in practice. 

It is the goal of the present paper to give a practical algorithm that 
actually finds the solutions if the coefficients and degrees of F are not pro¬ 
hibitively large. We believe that a pleasant feature of our algorithm is, that 
we do not use Puiseux series (only truncated power series) and that we work 
entirely over Q. 

For the sake of completeness we formulate Runge’s Theorem in a more 
algebraic geometric language and which also works in number fields. We 
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learnt this formulation from an informal note by Yu.Bilu. Let C be a smooth, 
connected algebraic curve defined over a number field K. Let S' be a finite 
set of places of K, including the infinite ones and let Os be the ring of 
S-integers. Fix a function / £ C{K). A point P £ C(K) will be called 
S-integral with respect to / if f(P) £ Os- The Galois group Gal(K/K) 
acts on the set of poles of /. Denote by E the set of orbits under this action. 

Theorem 1.2 (Runge) . Assume that |E] > |S|. Then the S-integral 
points of C are effectively bounded. 

For example, when K = Q, / is the x-coordinate and S consists of the 
place at oo, we see that we must have |E| > 1, i.e. the set of poles of x 
must exist of at least two Galois orbits. This is precisely the factorisation 
condition discussed earlier. Instead of taking one function / we could have 
taken a finite set of functions f\..... ft. In our case over Q we would take 
fi =x,f 2 = y- 

Using ideas of Sprindzuk, Bombieri [Bo] found an interesting extension 
of Runge’s theorem. Let s be a positive integer. A point P £ C(K) is called 
s-integral if \f(P)\ v > 1 for at most s places v of K. 

Theorem 1.3 (Bombieri, Sprindzuk) Assume that |E| > s. Then the 
s-integral points of C are effectively bounded. 

2 Preparations 

To start with, we assume that the Newton polygon of F has a slope E which 
is neither vertical nor horizontal (we call this a tilted slope). The remaining 
case, when the Newton polygon is a rectangle, will be dealt with at the end 
of this section. 

When there is only one slope, we assume that the associated polynomial 
Fe is a product of two relatively prime polynomials. This is called the 
Runge assumption. When there are only two slopes and one is horizontal 
and the other E we interchange x, y. We then get a new polynomial F whose 
Newton polygon has a vertical slope and a tilted slope E. After having made 
this change, if necessary, we are now in the position that Fe factors into 
two relatively prime polynomials whose degrees in y are strictly less than 
deg y (F). 

Suppose that the points on the slope E satisfy ax + by = w where a, b 
are relatively prime integers and w is some integer. We define the weight 
of a monomial x m y n by am, + bn. The weight of a polynomial P is the 
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maximum of the weights of the monomials occurring in P. Notation: w{P). 
In particular we have that w = w(F). 

Let us denote the factorisation of Fe by Fe = GqHq. We now give 
an algorithm to solve F(x,y) = 0. First we consider the real points of 
F(x,y ) = 0. Let x% be the largest positive zero of the discr J/ (F) and x\ 
the smallest. Take X 2 = x\ to be specified later if there are no real zeros. 
For any x > X 2 the equation F(x, y) = 0 has a fixed number, say r, of real 
solutions that we denote by y\, y2, ■ ■ ■, y r - We call the functions yi{x) the 
positive real branches of the curve F = 0. Note that for any branch yi(x) 
there is a real zero a* of Fo(l,a) such that yi{x)/x b ^ a —> a% as x —> 00 . As 
a consequence of the Runge assumption we can find a polynomial P{(x, y) 
with integer coefficients, not divisible by F, such that Pi(x,yi(x)) —> 0 as 
x —> 00 . The construction of Pi will be carried out in the next section. Here 
we conclude the algorithm. Choose a positive parameter r. Let x + {i ) be 
the largest positive zero of the resultant of F(x,y ) and Pi(x,y ) + r with 
respect to y. Let x~(i ) be the largest positive zero of the resultant of F and 
Pi(x,y ) — r. Let X % be the maximum of X2, x + (i), x~ (i) where we ignore 
X 2 if it has not been defined yet. Then we know that for all x > X,; we 
have \Pi(x,yi(x))\ < r. Suppose we have an integer point (x,y) on the i- 
th branch with x > X{. Then P%(x, y) assumes an integral value a with 
|a| < t. We find all such (x,y) simply by solving the simultaneous systems 
F(x,y) = 0 ,Pi(x,y) = a for all integers a with |a| < r. 

We carry out these steps for each positive real branch of F = 0. After 
that we have found all integer solutions (x,y) with x > max,; X t . Next we 
should consider the case x < x\. For any such x the equation F(x, y) = 0 has 
a fixed number r' of real solutions y%{x) which we call the negative branches. 
For each such i we construct a function Pi(x,y ) and proceed as above. 

In practice the number of distinct Pj is smaller than the number of actual 
real branches because one polynomial may vanish on several asymptotic 
branches of the curve. 

Finally, we promised to give an algorithm in the case when the Newton 
polygon has no tilted slopes, i.e. it is a rectangle. In that case F contains a 
term ax m y n where m = deg X (F) and n = deg y (F). When F(x,y ) = 0 with 
x large, the value of y will be close to a zero of Fziy) = lim^oo x~ m F(x, y). 
We now simply use the above algorithm by taking Fiiy) as polynomial whose 
value tends to zero as we let (x, y) follow a branch of F = 0 with x —■»■ 00 . 
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3 Construction of the vanishing functions 


Let notation be as in the introduction. Suppose we want to compute the 
integral points on one of the positive real branches y/x b ^ a —> a of F = 0. In 
this section we construct a polynomial whose value on this branch tends to 
0 as we let x —> oo. 

Perform the following change of variables, x —> 1 /t b ,y —> g/t a . We 
obtain 


F(x,y) 


1 

t w ( F ) 




1 

t w ( F ) 


(/o + tfi + t 2 f 2 + • • •) 


where G Z[r/] and /o = Go(l,rj)Ho(l,r]). We now perform a Hensel lift of 
/o(r?) = Go(l,r])Ho(l,r]) to a factorisation in (Q[[i]])[y] of the form 


f(t, rj) = ( g 0 +git-\ - )(h 0 + hit H-) 


where go = Go (1, rj) , ho = Ho(l,rj), the degrees of g t , i > 0 are strictly less 
than the degree of go and the degrees of hi, i > 0 are strictly less than the 
degree of ho- 

Notice that g := go(y) + gi(y)t + • • • = 0 is an analytic curve which, for 
small t, contains a subset of the branches of F = 0. Suppose our particular 
branch is among this union of branches. We now construct a polynomial P 
which vanishes on the branches of g = 0 as x —» oo. Consider a polynomial 
P € Z[x,y] with unknown coefficients and such that deg P < deg y (F) and 
w(P) < N for some integer N to be specified. We rewrite P{\/t a ,y/t b ) = 
t~i v ( p )p(t, rj). We choose our coefficients such that p (mod g) = 0(t N+1 ). 
The number N is chosen in such a way that the number of coefficients of P 
exceeds the number of equations following from the constraint. 


Lemma 3.1 The vector space 

{P <E Q[x,y] | deg y P < deg y F, w(P ) < N} 

has Q-dimension at least N' 2 /(25 x 5 y ) + N/{28 x ) + N/(2S y ) + C\{F) if N < 
w(F) and dimension at least Ndeg y (F)/a + (^(T 1 ) if N > w(F). Here 
C \, 6 b depend only on F, not on N. 

As to the number of equations provided by p (mod g) = 0(t N+1 ), a priori 
we have Ndeg v (go) conditions. However, this number would in general be 
too large. Fortunately we have the following additional consideration. Let 
C be a primitive a-th root of unity. Replacing t by ft and 7/ by r)( b does 
not change P(l/t a ,r]/t b ). Hence p(t,r 7 ) changes by a factor f~ b . A similar 
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remark holds for g. Hence, after computation of p (mod 5 ), the only terms 
that occur transform with the same character under our substitution. So 
the actual number of constraints is at most Ndeg v (go) / a + 1 . 

Because deg, ? (go) < deg y (F) we have that Ndeg y (F)/a + C 2 exceeds 
Ndeg v (go)/a for sufficiently large N. Hence there exists a polynomial P of 
the required type. 

4 Example 1 

Consider the equation 

F(x, y ) := y 6 — 2 y 5 — 4y 2 x 4 + 17yx 2 + 4x — 18 = 0. 

The highest degree part is given by y 6 —4y 2 x 4 . We now replace x —> 1 /i, y —► 
rj/t to get 

f(t, rj) := 4f 5 - 18t 6 + 17 t 3 g - 4r / 2 - 2 trf + rf = 0. 

We shall be interested in Hensel lifts of the factorisation 

if — 4 if = (if — 2 )(rf + 2 )rf 

up to order 4. We get 

= yiy2ff3 

where 

51 = if — 2 — gt — t 2 /2 + 15gt 3 /8 + 0(t 4 ) 

52 = if + 2 — ijt — t 2 /2 + 19i]t 3 /8 + O^ 4 ) 

53 = if - 17yt 3 /4 + 0(t A ). 

First we determine a polynomial that vanishes on the branches given by 
51 = 0. Let 

P(x, y) = (ax + by)(y 2 — 2x 2 ) + py 2 + qxy + rx 2 + kx + ly + m 
where a, b,p, q , r, k, l, m are numbers to be determined. Define 

p(t, rj) = t 3 P(l/t,p/t). 

Then, 

p (mod 51 ) = (26 + 2p + r + ay + qp)t 

+(a/2 + k + 3bi]/2 + Ip + pp)t 2 
+(—136/4 + m + p/2 — 15ap/8 )t 3 + 0(t A ). 
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This remainder vanishes up to order 4 if and only if 

2b + 2p + r = 0 a + q = 0 

a/2 + k = 0 3b/2 + l+p = 0 

— 136/4 + m + p/2 = 0 -15o/8 = 0. 

One solution is given by 

a = 0, 6 = 2,p = —3, q = 0, r = 2, k = 0, l = 0, m = 8 . 

So we get the desired function 

Pi = 2 y(y 2 - 2x 2 ) + 2x 2 - 3 y 2 + 8 . 

The integer points lying on g 2 = 0 do not correspond to any real branches 
extending to infinity because y 2 + 2x 2 has no real factors. Finally, the 
function y 2 vanishes on the branches given by = 0 as t —> 0 . 

Let us now solve the equation F(x,y ) = 0. First of all the real roots of 
the discriminant of F with respect to y lie in between 1 and 1.25. First we 
deal with the branches given by g\ = 0. The real zeros of the resultant of F 
and Pi + 1 with respect to y he between —4.1 and 3.2. The real zeros of the 
resultant of F and Pi — 1 are between —3.8 and 3. Hence, for the solutions 
on the branches given by g\ = 0 we have that |Pi(x, y)\ < 1 whenever x > 4 
or x < —5. In other words, we have P±(x,y ) = 0 for such points. The 
y-resultant of F and P± has no integer zeros x. 

As we said we can ignore the factor <72 because it does not correspond 
to ant branches. Finally we consider the branches given by 53 = 0. The 
resultant of y 2 + 1 and F has no real zeros, the resultant of y 2 — 1 and F 
has its real zeros between 1 and 2 . So y 2 < 1 for all integer points on the 
branches given by 53 = 0. In other words, y 2 = 0 for such points. Since 
y = 0 and F(x, y) = 0 imply 4x — 18 = 0, there are no integer solutions. 

We are left to check the remaining values of x between 4 and —5. So we 
check F(k, y) = 0 for integer solution for k = —4, —3,..., 2, 3. It turns out 
that there are no integer solutions. 

5 Example 2 

Consider the equation 

F(x, y) := y 4 + 2 y 3 — 9 x 2 y 2 + 2 xy — 15x — 7 = 0. 


7 



The highest degree part is given by y 4 — 9 x 2 y 2 = (y — 3 x)(y + 3 x)y 2 . First 
we define 

f{t,rj) = t 4 F(l/t,y/t). 

Then factor / up to order 4 as 

f(t,rj)= 515-253 

where 

5 ! = ?7-3 + t-t 2 /18- 13t 3 /54 + 0(t 4 ) 
g 2 = v + 3 + t + 5t 2 /18-13t 3 /54 + 0{t 4 ) 

53 = if ~ 2??t 2 /9 + 5t 3 /3 + 0(t 4 ). 

Clearly P ± := y — 3x + 1 vanishes on the branch given by 51 =0 and 
P '2 := y + 3x + 1 vanishes on the branch given by 52 = 0. A straightforward 
computation as in the previous section shows that P 3 := 2y 3 + 15y 2 vanishes 
on the branches given by 53 = 0 . 

We summarise the resultant computations in the following table. 


Resultant^ 

^min 

2-max 

F,Fy 

-1.03 

0.78 

F.P+l 

none 

none 

F,Pi- 1 

-1.46 

1.52 

F, P 2 + 1 

-1.15 

1.10 

PP 2 -1 

- 0.86 

0.75 

F,P3 + 1 

-2.18 

2.12 

P,P s -3 

-7.83 

2.10 


For the branches given by 51 = 0 and 52 = 0 we see that |Pi(x,y)| < 1 
and \P- 2 (x,y)\ < 1 for all integer points with \x\ > 2. In other words, 
Pi(x,y) = P- 2 (x,y) = 0. The system Pi = F = 0 does not give integer 
solutions, neither does P 2 = F = 0. 

For the branches given by 53 = 0 we see that —1 < P$(x, y) < 3 for 
all x with x < — 8 and x > 3. Hence P^(x,y) = 0,1,2. Combining these 
possibilities with F = 0 again does not yield any integral solutions. 

It remains to check all points with —7 < x < 2. When x = —1 we get 
the solutions y = —4, —1,1, 2. 

6 Example 3 

Consider the equation 

F(x, y) := (y 2 — x 3 )(y 2 — 2x 3 ) + 2x 5 — 9 xy — 3 = 0. 
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This is an example which satisfies the Runge condition with respect to the 
weight 2m + 3 n for a monomial x m y n . Introduce 

f(t,r ] )=t 12 F(l/t 2 ,r ] /t 3 ). 

We can Hensel lift this to a factorisation 

= am 

where 

gi = rf - 1 - 2 f 2 - 4t 4 - 16t 6 + 0(t 7 ) 
g . 2 = , ? 2 -2 + 2t 2 + 4t 4 + 16t 6 + 0(f 7 ). 

We look for a function of the form 

P := ay 2 + by + cyx + px 3 + qx 2 + rx + s 
which vanishes on the branches given by <71 = 0. Define 

P(t,v) ■= t 6 P(l/t 2 , rj/t 3 ). 

Then 

p (mod gi) = (a + p) + crjt + ( 2 a + q)t 2 + br/t 3 
+(4 a + r)t 4 + (16a + s)t 6 + 0(t 7 ). 

Note, by the way, that only terms t 2k ,pt 2l+1 occur. To have this remainder 
vanish up to order 7 we can take 

a = —1, b = 0, c = 0,p = 1, q = 2, r = 4, s = 16 

and we get the function 

P\ = x 3 — y 2 + 2x 2 + 4x + 16. 

Similarly we have 

p (mod < 72 ) = 2 a+p + crjt + (—2 a + q)t 2 + br/t 3 

((—4a + r)f 4 + (—16a + s)t 6 . 

To have this remainder vanish up to order 7 we can take 

a = 1, b = 0, c = 0,p = —2, q = 2, r = 4, s = 16 
and we get the function 

P -2 = y 2 — 2 x 3 + 2 x 2 + 4x + 16. 

We now make our table of ranges of real zeros for the various resultants. 
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Result ant y 

2 -min 

2 -max 

F,Fy 

-1.50 

8.23 

F, Pi + 5 

15.0 

31.7 

F, P 1 - 0.5 

none 

none 

F, P-2 + 5 

13.8 

36.1 

f,p 2 - 1 

none 

none 


So, when x < — 2 or x > 37 we have either —5 < P\{x,y) < 0.5 or —5 < 

Piix, y) < 1. First we solve P\(x, y) + k = 0, F(x, y) = 0 for k = 0,1, 2 , 3,4. 

There are no integer solutions. Then we solve P 2 (x,y) + k = 0 ,F(x,y) = 0 

for k = 0,1,2,3,4. Again there are no solutions. Finally we solve F(k,y) = 0 

for — 1 < k < 36. We find the solution x = 2, y = 3. 
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